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ABSTRACT 

An  analytical  investigation  was  made  of  combined  forced 
and  free  turbulent  convection  in  a  circular  vertical  tube  with 
volume  heat  sources  and  constant  wall  heat  addition. 

The  IBM  1620  digital  computers  located  at  The  George  Wash- 
ington University  and  also  at  the  United  States  Naval  Academy 
were  used  to  solve  the  basic  equations  in  which  the  parameters 
are:   Prandtl  number  Pr;  Rayleigh  number  Ra;  friction  Reynolds 
number  Re*;  and  the  volume  heat  source  parameter  P.   For  fixed 
values  of  these  parameters  the  solution  gave  the  fully  developed 
velocity  profile,  the  temperature  profile,  and  the  pressure  drop. 
The  program  also  solved  for  the  mixed-mean-to-wall  temperature 
difference,  Nusselt  number,  and  Reynolds  number.   The  following 
values  of  these  parameters  were  investigated:   Pr  =  1,  10,  100; 
Ra  =  0,  3^,  5^,  10^;  Re*  =  103,  10^,  105;  and  F  =  .1,  .5,  1,  10. 

The  results  of  laminar  flow  problems  were  investigated  to 
check  the  validity  of  the  program  and  they  compared  favorably 
with  known  previous  results.   In  the  case  of  pure  forced  convec- 
tion, laminar  flow  with  a  negligible  volume  heat  source,  the  ve- 
locity profile  and  the  Nusselt  number  checked  exactly  with  Dr. 
Ojalvo's  results  for  the  same  problem  with  no  volume  heat  sources 


The  results  for  turbulent  heat  transfer  in  combined  forced 
and  free  convection  with  volume  heat  sources  showed  that  an  in- 
crease in  Prandtl  number  had  a  smaller  effect  than  an  increase 
in  Re*  on  increasing  Nusse.lt  number.   The  volume  heat  source 
parameter  P  had  negligible  effect  on  Nu„  Rayleigh  number  had 
pa. effect  on  Nu.   The  pressure-drop  parameter-  C  increased  approx- 
imately an  order  of  magnitude  as  Re*  increased  an  order  of  mag- 
nitude.  Prandtl  number  Pr,  Ra,  and  P  had  no  effect  on  C  for 
Ra  <  625.  When  Ra  >  625  and  F  >  1,  increasing  Pr  decreased  C 
and  increasing  Ra  increased  C;  and  for  Ra  >  625  and  F  <  1,  in- 
creasing Pr  increased  C  and  increasing  Ra  decreased  C.   Decreas- 
ing P  for  Ra  >  625  decreased  C.   Increasing  Pr  had  less  effect 
than  increasing  Re*  on  decreasing  the  mixed--mean~to~wall  temper- 
ature rdiff erence  cj>rn  .   Increasing  Ra  had  negligible  effect  on 
(f»rn  ,  Decreasing  P  increased  c|)rn  . 

The  velocity  in  the  center  of  the  tube  lowered  and  became 
less  positive  as  Pr,  Re*,  or  volume  heat  sources  were  Increased. 
When  P  <  1,  increasing  Ra  flattened  the  velocity  profile. 

Increasing  Re*  or  Pr  flattened  the  temperature  difference 
profile  (j)  .   Increasing  Ra  had  negligible  effect  on   <j)   .   De- 
creasing P  increased  (b   and  when  F  =  1,  the  temperature  differ- 
ence profile  passed  through  zero. 
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CHAPTER  I 
INTRODUCTION 

The  first  theoretical  investigations?  Into  combined  forced 
and  free  laminar  flow  in  a  vertical  tube  were  conducted  by 
Ostroumov-^"  and  Hallmark ,   in  addition  to  predicting  velocity 
and  temperature  profiles,  Hallman  extended  the  analysis  to  in- 
clude calculation  of  Nusselt  numbers  and  pressure  drops.   He 
considered  cases  in  which  volume  heat  sources  were  either  pres- 
ent or  not  present.   He  also  conducted  an  experimental  investi- 
gation^ ..of  the  problem  for  laminar  flow. 

The  first  theoretical  investigation  into  combined  forced 
and  free  turbulent  convection  in  a  vertical  tube  was  conducted 
by  Ojalvo  and  Grosh".   The  analysis  did  not  consider  cases  with 
volume  heat  sources  present.   All  cases  of  laminar  heat  trans- 
fer, pure  forced  convection  or  combined  forced  and  free  convec- 
tion, checked  exactly  with  known  results.   Cases  of  pure  forced 
convection  turbulent  flow  were  net  in  agreement  with  known  re- 
sults apparently  because  too  high  a  value  of  eddy   viscosity  was 
used  in  the  buffer  zone  and  in  part  of  the  turbulent  core.   The 


^'Superscript  numbers  refer  to  similarly  numbered  "references  in 
the  Bibliography. 


results  of  combined  forced  and  free  turbulent  convection  were 
presented. 

Sackett^-3  presented  data  for  combined  forced  and  free  tur- 
bulent convection  in  a  vertical  tube  with  a  new  relationship  for 
the  eddy  diffusivity  of  momentum.   The  results  for  turbulent 
flow  cases  indicated  some  improvement  but  did  not  compare  favor- 
ably with  Dr.  Ojalvo's  results.   The  most  likely  cause  was  attri- 
buted to  the  reworked  computer  program. 

Jackson^  conducted  an  investigation  to  improve  the  relation- 
ship for  the  eddy  diffusivity  of  momentum.   The  results  of  his 
study  are  used  in  the  present  study. 

The  present  analysis  investigates  the  transfer  of  uniform 
thermal  energy  from  the  wall  of  a  round  tube  and  uniform  ther- 
mal energy  from  volume  heat  sources  to  a  fluid  flowing  vertically 
upwards  in  the  tube.   Both  forced  and  free  convection  exist  in 
the  steady  turbulent  flow.   Only  the  case  of  fully  developed  flow 
is  considered  and  compared  with  known  results,, 

In  reworking  the  original  computer  program,  an  error  was 
found  that  would  affect  Dr.  Ojalvo's  results  for  combined  forced 
and  free  turbulent  convection.   This  error  did  not  occur  in 
Sackett's  program.,  thus  the  apparent  cause  for  the  differences 
in  their  results. 


The  analysis  has  practical  applications  in  the  fields  of 
nuclear  reactors  and  heat  exchangers.   Its  solution  will  show 
the  effect  of  uniform  volume  heat  sources  on  combined  forced 
and  free  turbulent  convection  problems. 


CHAPTER  II 
ANALYSIS 

The  problem  to  be  analyzed  is  combined  forced  and  free 
turbulent  convection  with  uniform  volume  heat  sources  in  a 
vertical  circular  tube  whose  axis  is  parallel  to  the  direction 
of  the  body  force.   There  is  to  be  a  net  through- flow. 

In  addition  to  the  description  of  the  problem  given  above, 
the  following  assumptions  are  made: 

1.  Axial  symmetry  exists  for  the  momentum  and  heat 
transfer. 

2.  All  fluid  properties  are  constant,  except  density 
in  the  expression  for  body  force.  A  mean  density 
Pf*  *  is  used  for  all  other  density  terms. 

3.  Viscous  dissipation  and  axial  heat  conduction  are 
negligible  compared  with  the  heat  conduction  in 
the  radial  direction. 

4.  There  is  uniform  wall  heat  flux, 

5.  The  velocity  and  temperature  profiles  are  fully  de- 
veloped.  There  are  no  radial  or  angular  velocity 
components. 

6.  There  is  single-phase  flow. 
•The  Nomenclature  is  given  in  Appendix  A. 


7.  Flov/  is  steady j  turbulent  and  incompressible. 

8.  The  eddy  diffusivity  of  momentum  is  given  by  Jack- 
son's equations.   (See  Appendix  B.) 


The  basic  equations  employed  in  this  analysis  are  the 
continuity,,  Navier-Stokes,  and  energy  equations  in  cylindrical 
coordinates^ K      An  equation  of  state  is  also  used, 

On  the  basis  of  the  above  assumptions  and  description  of 
the  problem,  the  continuity  equation 

It    +   V.  (pvj    =    0 

reduces  to 


&  =  °   > 


Wti'f  =  hMr[k+^»]& 


a) 


the  Navier-Stokes    (Momentum)    equations 

^  4|   =  B  ^  -  Vp  iy^V2  V  f  y/<  V  (V  •  V  ) 

reduces   to 

a*   ^Sb  r   arV   Life  $c  J  *»*' 

and  the  energy  equation 
reduces  to 


3f-/  .       (4) 
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The  equation  of  state  is  developed  from  an  expansion  of  o     in 

a  Taylor  series  in  t  about  the  reference  temperature  tw,  and 
8 


is: 


(5) 


?    =   ?w[l-0(Hw). 

J^SX? lpPfleflfc   of  Equ_a tion s 

Equation  (l)  indicates  that 

u  =  u  (r),  (6) 

and  equations  (3)  indicate  that 

P  =  P  (x).  (7) 

If  equations  (5),  (6),  and  (7)  are  utilised,  equation  (2)  can 

>• 
be  written  as 


1 d 

r 


(8) 


The  boundary  condition  of  uniform  wall  heat  flux  plus 
the  assumptions  of  constant  specific  heat,  uniform  volume  heat 
sources,  and  a  fully  developed  temperature  profile  require  that 


R 


(a  constant) . 


(9) 


Equation  (9)  is  developed  in  Appendix  C.   A  new  variable  0 
is  introduced  and  defined  as 


&  =  0Ci-)  =  t(x,i-)-t(x,-£)  =  j-tw 


(10) 


In  terms  of  equation  (10),  equations  (4)  and  (8)  become, 


respectively, 


Ucr  u  R  -  f 


I 


r    lr 


k+tm*'»]$ 


(id 


and 


*■  i'r 


Uit 


M   t     flH 


ftWtfwPf  e<« 


#  fx) 


r  \v/ 


9c 


(12) 


Although  the  density  P  t      in  equation  (12)  Is  actually  a 
function  of  x,  by  considering  the  solutions  of  equations  (11 ) 
and  (12)  to  obtain  u  and  G  at  a  fixed  value  of  x,  we  can  con- 
sider Pw  constant".   Our  assumption  of  fully  developed  flow 
and  heat  transfer  is  consistant  with  these  results.   Each  side 
of  equation  (12)  may  be  set  equal  to  some  constant  since  we 
have  separated  the  variables^.   Thus, 


¥  t  p  - 


9 
3c 


3llA  rnjU  C 

Da3c 


(  a   const  a  n  t ) 


(13) 


and 


r  dr 


r 


^^fffi)^^l^--'triI--(i4) 


u   3c 

The  pressure-drop  parameter  C  in  these  equations  takes  on  the 
value  of  unity  for  the  special  case  of  pure  forced-convection 
laminar  low,  as  described  by  Hallman^. 


DO 


In  order  to  solve  the  set  of  problems  described.,  equations 
(11)  and  (14)  are  nondimensionalized  by  using  the  following 
dimensionless  terms: 


radius 


temperature  difference   (l6) 


velocity 


7)  e  £r/D 
$  =  lbk6/i"0? 
U  =  u/U  m 

Rq=fWfw?CPfli?M6/»k      Rayleigh  number 
F   =  fm^M    T    /   °\,  volume  heat   source 

Equation  (19)  is  discussed  in  Appendix  P. 

Thus j  equation  (ll)  becomes 


J- 4(7) 
7)  d7|V 


1  + 


Lt\ 


m 


44 

d7] 


1-FV-  f- 


and  equation  (l4)  become; 


Ufi 


■  Cm 


(fa  d>      or 


In  accordance  with  assumption  (8),  we  will  us* 
.in  -  fl A     <> / 

6  —mm  W         ~~*  *i       ^ 


1 


w 


I  Nze*-f>(.0\H*/Pr) 


_6 


(15) 


(17) 
(18) 

(19) 


(20) 


(21) 


(22) 


as  given  by  Lykoudis  .  Equation  (20)  may  thus  be  written  as 


">1  (/71  (^ 


1  + 


trrt- 


M 


7/_ 


d7| 


j  =  ^F  [/-</" 


(23) 
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Bound  a  ry  C  o  nd  i  t  ions 
The  following  boundary  conditions  are  used: 

Ud)    -    0  (2") 

$  0)  =  0  (25) 

jiJtLifi)  =  0'(o)  -  0  (27) 

These  boundary  conditions  come  from  the  physical  problem. 
Equations  (24)  and  (25)  state  that  the  fluid  velocity  and  tem- 
perature at  the  wall  of  the  tube  (7)  -  l)  are  equal,  respec- 
tively ,    to  the  wall  velocity  and  wall  temperature.   Equations 
(26)  and  (27)  come  from  the  fact  that  heat  and  momentum  are 
not  transferred  across  the  center  line  [y\   =  Q)    of  the  tube, 
due  to  axial  symmetry,  resulting  in  a  zero  slope  for  the  temper 
ature  and  velocity  profiles  at  the  center  line7. 

Equation  (21),  the  momentum  equation,  and  equation  (23) * 
the  energy  equation,  are  to  be  solved  by  a  1620  IBM  digital 
computer  for  U  and  (f)  as  functions  of  7)  .   All  boundary  con- 
ditions are  to  be  satisfied.   The  value  of  the  pressure-drop 
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parameter  C  may  also  be  obtained  if  we  use  the  following  form 
of  the  continuity  equation: 

Urn  =e/U"7|d"H=l    •  (28) 

Three  dimensionless  parameters  in  our  equations  are  the 
Prandtl  number  Pr^  the  Rayleigh  number  Ra>  and  the  volume  heat 
source  parameter  Fa      The  Prandtl  number  determines  the  rela- 
tionship between  the  velocity  and  temperature  distributions. 
For  the  special  case  of  Pr  a  1,  both  profiles  are  the  same. 
Rayleigh  number  Ra  z   Gr  Pr  where  the  Grashof  number  Gr  repre- 
sents a  ratio  of  buoyant  forces  to  viscous  forces^2.   Thus., 
free-convection  predominates  for  large  Rayleigh  numbers  and 
for  the  special  case  of  Ra  =  0.,  pure  forced-convection  is  rep- 
resented.  The  significance  of  variations  in  the  magnitude  of 
the  volume  heat  source  parameter  P  is  discussed  in  Appendix  P. 

Empirical  equations  for  -^7-   are  given  in  Appendix  B. 
They  are:  g  _ 

^z.ERe'V'-'n)   -|  forj^$y|<|-j^(30) 

^=    IffeVfOO  _|  for  0s<^l,      (3D 

where  a e   r    0  U^-^y  friction  Reynolds   number      /-^2) 

and  W^e    V?vv  9c/fvv     /    friction  velocity  (33) 
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Thus,  a  fourth  parameter,  the  friction  Reynolds  number  Re*, 
is  introduced.   In  this  analysis,  when  Re*  =  0,  the  ratio  — J 
is  set  equal  to  zero,  thus  assisting  in  checking  the  validity 
of  the  program  by  comparison  with  laminar  flow  problem  results. 

Exte2ision  of  Results 

The  solution  of  equations  (21),  (23),  and  (28)  will  yield 
U  and  (f)  as  functions  of  7}  for  a  given  set  of  values  of  the 
parameters:   Ra,  Pr,  Re*,  and  P.   The  pressure-drop  parameter 
C  will  also  be  obtained. 

These  results  may  be  extended  by  calculating  the  following 
useful  quantities:   the  dimensionless  mixed-mean-to-wall  tem- 
perature difference  <$m    ;    the  Nusselt  number  Nu;  and  the  Rey- 
nolds number  Re. 

Hallman3  gives  equations  for  the  first  two  of  these  quan- 
tities.  They  are: 

and 

Nu    E    -£-    (}~F)    ;  (35) 

see  Appendix  D. 

An  expression  for  the  Reynolds  number  is  derived  In 
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Appendix  E.   The  result  is 

|?e  5         ^(|e  >  ,  (36) 

where  the  -j-  sign  is  for  upward  flow  and  the  -  sign  is  for 
net  downward  flow. 
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CHAPTER  III 
METHOD  OF  SOLUTION 

The  method  of  solution  is  based  on  the  use  of  a  digital 
computer.   This  method  will  decrease  considerably  the  amount 
of  hand  calculations  required  in  the  solution  of  our  problems. 
The  possibility  of  errors  are  thus  reduced.   This  method  of 
solution  will  also  allow  us  to  investigate  many  different 
problems  and  compare  some  of  these  results  with  known  analy- 
tical and  experimental  data. 

The  independent  parameters  used  are  Ra,  Re*,  Pr,  and  P. 
The  Rayleigh  number  is  chosen  to  measure  the  extent  of  the 
free-convective  effect  on  the  flovj.   The  friction  Reynolds 
number  determines  the  value  of  -7^-  as  a  function  of  7]  . 
The  Prandtl  number  determines  the  value  of  <T"  .   The  volume 
heat  source  parameter  is  chosen  to  describe  the  thermal  ener- 
gy convected  downstream  to  the  heat  generated  in  the  fluid. 
V/ith  these  input  quantities  fixed,  the  method  of  solution  pro- 
ceeds as  follows: 

1.   Ui  =  2  (1  -7J2)  (37) 

is  assumed  as  an  initial  guess  for  the  velocity  profile.   This 
equation  satisfies  the  continuity  equation   /  U^ld^l     -  _L 
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and  the  boundary  conditions  U(l)  ~  0,  and  U' (0)   -  0.   It  is 
the  limiting  case  of  the  parabolic  profile  in  pure  forced- 
convection  laminar  flow. 

2.  Equation  (23)  is  integrated  with  the  aid  of  boundary  .«' 

i 

condition    (jj    ($)=.   0  to  obtain 

a/~    jJLJihliJichi   _    2_2i .        r38) 

V"        7\[|  +  (rPri^J  [|+<rR-|J 

/^  —J 

3.  Equation   (38)    is  integrated  using  boundary  condition 
(j)  (l)   r  0  to  obtain 

<b     =    p  ty'd1\  (39) 

as  a  function  of  ~Y\    . 

4.  Equation  (21)  is  integrated  using  boundary  condition 
U1  (0)  =  0  to  obtain 

"       &/  =  _    li2LiVU23 f_2LC .    (40  j 

The  pressure-drop  parameter  C  must  have  an  output  numer- 
ical value  for  the  computer.  For  the  limiting  case  of  pure 
forced-convection  laminar  flow  C  =  1. 

5.  Equation  (ho)   is  integrated  subject  to  boundary  condi- 
tion U(l)  =  0  to  obtain 

Ut   --pvU-r\  (41) 

as  a  function  of  ~Y)     . 

6.  The  continuity  equation  (28)  is  cheeked  to  see  if 
equation  (;-!-l)  satisfies  it.  If  U2  dees  not  satisfy  this  equation, 
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the  value  of  C  is  changed  until  it  is  satisfied. 

7.  The  function  U2  ^s   checked  against  U'x  by  comparing 
U'gC1)  with  U'^fl)  to  see  if  they  are  within  0.1$  of  each  other 

8.  Equations 


m 


l/u 


'! 


(I-F) 


m 


*\2. 


and  Re  z 


t  ( j?e*) 


8  C  +  /?a  ^/f  p 

are  used  to  calculate  d)rVl  ,  Nu  and  Re,  respectively ,  if  the 
check  in  Step  7  is  satisfactory.   The  results  are  then  printed 
out  as  given  in  Step  10  below. 

9.   If  the  check  in  Step  7  is  not  satisfactory,  U2  is  used 
as  an  initial  assumption,  replacing  U'x  in  going  through  the 
procedure  again,  starting  at  Step  2. 
10.   Print  out 


Pr ob 1 em  Numb  e r 


Ra 
1 


Re 

u2 

v 

0 


4> 


^ 


Title 

Pr 
Ul 

V 

0 

Nu 


F 

(T 

v 

0 

Re 
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Up  and  Ujl  are  both  printed  out  to  ensure  that  the  criterion 
for  checking  them  in  Step  7  is  valid. 

These  steps  may  be  shown  in  a  simplified  flow  diagram  for 
the  calculation,  Figure  1 „ 
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CHAPTER  IV 


COMPARISON  OP  RESULTS 


Table  1  lists  the  following  calculated  results:  0,   (j)  m, 
Re,  and  Nu  for  all  problems  solved  by  the  1620  digital  computer, 


TABLE  1 


RESULTS  FOR  ALL  COMPLETED  PROBLEMS 


Input  Parameters 
Problem 

No.  Ra  Re* 


Calculated  Results 


Pr 


(P 


m 


Re 


Nu 


403B   0    0 


A.  Laminar  flow,  pure  forced  convection 


103 


1.0  -915.85 


4.3 


o 


B.  Laminar  flow,  combined  forced  and  free 


40  2D 

81 

0 



10 

2.4 

-9.62 

-- 

3.7 

C.  Tu 

rbul en 

tflow, 

pure  forced 

convection 

379 

0 

103 

1 

.5 

7.3 

.036 

17000 

55.4 

380 

0 

10^ 

1 

1.0 

7.3 

-.0034 

17000 

0.0 

383 

0 

103 

1 

10.0 

7.3 

-.71 

17000 

50.3 

385 

0 

103 

10 

.5 

7.3 

.01 

17000 

196.0 

386 

0 

103 

10 

1.0 

7.3 

-.00053 

17000 

0.0 

389 

0 

103 

10 

10.0 

7.3 

-.194 

17000 

185.0 

391 

0 

103 

100 

.5 

7.3 

.002 

17000 

983.0 

392 

0 

103 

100 

1.0 

7.3 

«. 000064 

17000 

0.0 

395 

0 

103 

1. 

100 

10.0 

7.3 

-.038 

17000 

952.0 

397 

0 

102* 

1, 

1 

.5 

53.9 

.005 

232000 

4oo.o 

398 

0 

102* 

] 

1.0 

53.9 

-.00017 

232000 

0.0 

400B 

0 

300 

1 

10.0 

2.8 

-2.07 

4030 

17.4 

4oi 

0 

10* 

1, 

1 

10.0 

53.9 

-.093 

232000 

385.0 

403 

0 

10^ 

10^ 

1. 

10 

.5 

53.9 

.0011 

232000 

1750.0 

404 

0 

10 

1.0 

53.9 

-.000019 

232000 

0.0 

407 

0 

10 

10.0 

53.9 

-.02 

232000 

1720.0 

409 

0 

104 

100 

.5 

53.9 

.00025 

232000 

7900.0 
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Table  1  ( C  on  t  inu  eel ) 

Inpu t  Parame  t  e r* s 
Problem 

No.  Ra  Re*   Pr   P 


Calculated  Results 
(j)m       Re       Nu 


410 
413 
415 
416 

419 
421 
422 
425 

427 
428 
429B 


C. 
0 
0 

0 
0 
0 
0 
0 
0 
0 
0 
0 


Turbulent  flow., 
100 


10 
101* 

ior> 
105 
105 
105 
io5 
105 
105 

10? 
10^ 


100 

1 
1 
1 

10 

10 

10 

100 

100 

1 


1.0 
10.0 

.5 
1.0 

10,0 

.5 
1.0 

10.0 

.5 
1.0 
103 


>ure  forced  convecti 
53.9   -.000002 
53.9  -.0046 
430.0   .00063 
430.0   -.000021 
430.0   -.012 
430.0   .00017 
430.0   -.0000034 
430.0   -.0031 
430.0   .0001 
430.0  -.0000014 
53.9^10.34 


on  (Cont 

232000 

232000 

2900000 

2900000 

2900000 

2900000 

2900000 

2900000 

2900000 

2900000 

232000 


inued) 

0.0 

785O.O 

3150.0 

0.0 

3040.0 

11700.0 

0.0 

11500.0 

18600.0 

0.0 

386.0 


55 
56 

57 
60 
61 
62 

63 
66 
61 
68 

69 
72 
73 
74 
75 
78 

79 
80 
81 
84 

85 

86 

87 

90 


I).  Turbulent  flow, 
1    .1 


81  103 

81  103 

81  103 

81  103 

81  103 

81  103 

81  103 

81  103 

81  103 

81  103 

81  103 

81  103 

81  10? 

81  107 

81  10? 

81  107 

81  iof 

81  107 

81  107 

81  107 

81  10/ 

81  10| 

81  107 

81  104 


1  .5 

1  1.0 

1  10.0 

10  .1 

10  .5 

10  1.0 

10  10.0 

100  .1 

100  .5 

100  1.0 

100  10.0 

1  .1 

1  .5 

1  1.0 

1  10.0 

10  .1 

10  .5 

10  1.0 

10  10.0 

100  .1 

100  . 5 

100  1.6 

100  10.0 


combi 
5.8 

7.1 

7.3 

7.5 

6.9 
7.3 

7.3 

7.4 

7.2 

7.3 

7.3 

7.3 

53.7 

53.9 

53.9 

53.9 

53.8 

53.9 
53.9 
53.9 
53.9 
53.9 
53.9 
53.9 


ned  forced 

.067 

.036 
-.0036 
-.718 

.018 

.01 
-.00054 
-.194 

.0037 
.002 
-.000064 

-.037 

.  0091 

.005 
-.00017 

-.093 

.002 

.001 
-.00002 
-.021 

.00045 

.00025 
-.000002 
-.0045 


and  free 

1.6600 

17000 

17000 

17000 

I69OO 

17000 

17000 

17000 

17000 

17000 

17000 

17000 

232000 

232000 

232000 

232000 

232000 

232000 

232000 

232000 

232000 

232000 

232000 

23  2000 


53.1 

55.4 

0.0 

50.0 

191.0 

196.0 

0.0 

185.0 
970.0 
983.0 

0.0 

952.0 
394.0 

400.0 
0.0 

385.0 
1740.0 
1.750.0 

0.0 

1720.0 

7890 . 0 

7900 . 0 

0.0 

7850.0 
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Tabl 

e  1  ( 

Continued) 
Input  Parameters 

Calculated 

Results 

Problem 

No 

1.  Ra 

.  Re* 

Pr 

p 

c 

*m 

Re 

Nu 

D. 

Turbul 

ent  flow,  c 

ombined 

forced  and 

free  (Continued) 

91 

81 

105 

1 

.1 

427.0 

.0011 

2920000 

3120.0 

92 

81 

1 

.5 

430.0 

.00063 

2900000 

3150.0 

93 

81 

kP 

1 

1.0 

430,0 

-.000021 

2900000 

0.0 

96 

81 

icP 

1 

10 . 0 

430.0 

-.012 

2900000 

3040.0 

97 

81 

10| 

10J 

10 

.1 

427.0 

.  0003 

2920000 

11700.0 

98 

81 

10 

.5 

430,0 

.00017 

2900000 

11700.0 

99 

81 

105 

10 

1,0 

430 . 0 

-.0000034 

2900000 

0.0 

102 

.  81 

io5 

10 

10.0 

430.0 

-.0031 

2900000 

11500.0 

103 

81 

105 

100 

.1 

427,0 

.00019 

2920000 

18600.0 

104 

81 

105 

100 

.5 

430.0 

.0001 

2900000 

18600.0 

105 

81 

105 

100 

1.0 

430.0 

-.0000014 

2900000 

0.0 

108 

81 

105 

100 

10.0 

430.0 

-.0019 

2900000 

18400.0 

164 

625 

10^ 
10J 

1 

.5 

6.0 

.036 

16700 

55.4 

165 

625 

1 

1.0 

7.4 

-.0035 

17000 

0.0 

I69 

625 

10J 

10 

.1 

4.0 

.019 

16200 

191.0 

170 

625 

103 

10 

.5 

7.0 

.010 

17300 

196.0 

171 

625 

103 

10 

1.0 

7.3 

-.00056 

17000 

0.0 

174 

625 

103 
10j 

10 

10.0 

7.6 

-.194 

32000 

185.0 

175 

625 

100 

.1 

6.7 

.0037 

169OO 

970.0 

176 

625 

100 

.5 

7.2 

.002 

17100 

983.0 

177 

625 

10J 

100 

1.0 

7.3 

-.000064 

17000 

0.0 

180 

625 

10? 

10" 

10? 
10? 
10? 

10* 

1, 

100 

10.0 

7.4 

-..038 

18800 

952.0 

181 

625 

1 

.1 

52.2 

.0091 

231000 

394.0 

182 

625 

1 

.5 

53.7 

.005 

232000 

4oo.o 

183 

625 

1 

1.0 

53.9 

-.00017 

232000 

0.0 

I87 

625 

10 

.1 

53.5 

.002 

232000 

1740.0 

188 

625 

10 

.5 

53.8 

.0011 

232000 

1750.0 

189 

625 

lot 

10? 
10? 

10" 

1, 

10 

1.0 

53.9 

-.000019 

232000 

0.0 

192 

625 

10 

10.0 

53.9 

-.021 

233000 

1720.0 

193 

625 

100 

.1 

53.8 

.00046 

232000 

7890.0 

194 

625 

100 

.5 

53.9 

.00025 

232000 

7910.0 

195 

625 

10? 
105 

100 

1.0 

53.9 

-.000002 

232000 

0.0 

198 

625 

100 

10.0 

53.9 

-.0046 

232000 

7850.0 

199 

625 

1 

.1 

427.0 

.0011 

2920000 

3120.0 

200 

625 

105 

1 

.5 

430.0 

.00063 

2900000 

3150.0 

201 

625 

105 

1 

1.0 

430.0 

-.000021 

2900000 

0.0 

205 

625 

105 

10 

.1 

427.0 

.0003 

2920000 

11700.0 
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Table 

!  1  (c 

Jontinu 

ed) 

Input 

Parameters 

Cal 

culated  Re 

suits 

Probl 

.em 

No. 

Ra 

Re* 

Pr 

F 

c 

♦m 

Re 

Nu 

D.  Turbulent  fl 

ovj,,  combined 

forced  and 

free  (Cont 

inued) 

206 

625 

105 

10 

.5 

430.0 

.00017 

2910000 

11700.0 

207 

625 

105 

10 

1.0 

430.0 

-.0000034 

2910000 

0.0 

210 

625 

lo5 

10 

10.0 

430.0 

-.0031 

2910000 

11500.0 

211 

625 

105 
109 

100 

.1 

427.0 

.00019 

2920000 

18600.0 

212 

625 

100 

.5 

430.0 

.0001 

2910000 

18600.0 

213 

625 

lo5 

100 

1.0 

430.0 

-.0000014 

2910000 

0.0 

216 

62£ 
10? 
IO7 

1, 

lo5 

100 

10.0 

430.0 

-.0019 

2910000 

18300.0 

327 

103 

1 

1.0 

8.3 

"-.0035 

17400 

0.0 

330 

103 

0 

1 

10.0 

27.6 

-.73 

27600 

48.7 

332 

10J 

10 

.5 

1.6 

.01 

15700 

196.0 

333 

10 

1.0 

7.5 

-.00054 

17100 

0.0 

336 

1(F 

10 

10.0 

12.7 

-.194 

18700 

185.0 

338 

101( 

lo3 

100 

.5 

6.2 

.002 

I67OO 

984.0 

339 

107 

10? 

1, 

103 

100 

1.0 

7.4 

-.000072 

17000 

0.0 

342 

103 

100 

10.0 

8.4 

-.038 

17300 

948.0 

343 

10? 
10? 
10? 

'1 

.1 

27.4 

.0091 

223000 

393.0 

344- 

10? 

10h 

1 

.5 

51.1 

.005 

230000 

400.0 

345 

1 

1.0 

54.0 

-.0002 

231000 

0.0 

348 

10? 

10? 

10* 
4 

104 

104 
107 

10? 

10? 

10? 

10" 

1. 

10? 
10? 
10? 
10? 
10? 
10? 
10? 

104 
10- 

1 

10.0 

56.7 

-.094 

232000 

383.0 

349 

10 

.1 

47.8 

.002 

230000 

1740.0 

350 

10 

.5 

53.2 

.0011 

232000 

1750.0 

351 

10 

1.0 

53.9 

-.000019 

232000 

0.0 

354 

10 

10.0 

54.5 

-.021 

232000 

1720.0 

355 

100 

.1 

52.5 

.00045 

231000 

789O.O 

356 

100 

.5 

53.7 

.00025 

232000 

7910.0 

357 

100 

1.0 

53.9 

-.000002 

232000 

0.0 

360 

100 

10.5 

54.0 

-.0046 

23  2000 

7850.0 

361 

105 

1 

.1 

424.0 

.0011 

2920000 

3120.0 

362 

10* 

lo5 

1 

.5 

430.0 

.U0O63 

2900000 

3150.0 

363 

10? 

105 

1 

1.0 

430.0 

-.00006 

2900000 

0.0 

366 

10? 

10? 
10? 
10? 
10? 
10? 
10? 

10; 

10  y 

105 

1 

10.0 

430.0 

-.012 

2910000 

3050.0 

367 

105 

10 

.1 

426.0 

.00003 

2920000 

11700.0 

368 

lo5 

10 

.5 

430.0 

.00017 

2910000 

11700.0 

369 

105 

10 

1.0 

430.0 

-.0000034 

2910000 

0.0 

372 

105 

10 

10.0 

430.0 

-.0031 

2910000 

11500.0 

373 

105 

100 

.1 

427.0 

.00019 

2920000 

18600.0 

374 

10l 

100 

.5 

430.0 

.0001 

2910000 

18600.0 

375 

10l 

100 

1.0 

430.0 

-.0000014 

2910000 

0.0 

378 

105 

100 

10.0 

430.0 

-.0019 

2910000 

18400.0 
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Part  A  consists  of  the  problem  in  laminar  flow,  pure  forced 
convection  with  negligible  volume  heat  sources.   Part  B  con- 
sists of  a  problem  in  laminar  flow,  combined  forced  and  free 
convection  with  volume  heat  sources.   Part  C  consists  of  28 
problems  in  turbulent  flow,  pure  forced  convection  with  volume 
heat  sources.   Part  D  consists  of  100  problems  in  turbulent 
flow,  combined  forced  and  free  convection  with  volume  heat 
sources. 

The  IBM  solutions  to  all  problems  are  an  attachment  to 
the  present  analysis.   Copies  of  the  Fortran  program,  the  data 
cards,  and  the  solutions  to  all  problems  on  IBM  cards  are  in 
the  files  of  Dr.  Ojalvo, 

Discussion 

The  velocity  profile  and  the  Nusselt  number  of  the  laminar 

flow  problem  in  part  A,  Table  1,  checked  exactly  with  the  re- 
ft o 
suits  of  Ojalvo  and  Hallman3.   A  negligible  volume  heat  source, 

F=1000,  was  used  in  this  problem  to  check  the  validity  of  the 
program.   The  laminar  flow  problem,  combined  forced  and  free 
convection  with  volume  heat  sources  of  part  B,  Table  1,  com- 
pared favorably  with  Hallman0   Tt  also  compared  favorably  with 
Ojalvo  even  though  no  volume  heat  sources  were  present  in  his 
analysis. 

2k 


Problem  397  was  chosen  as  typical  of  the  problems  in  part 
C.   This  problem  was  compared  with  test  problem  429B  and  the 
universal  velocity  and  temperature  profiles  of  Eckert  and 
Drake2. 

The  values  of  U  and  (!)  m  in  problem  429B  with  negligible 

volume  heat  sources,  PrlOOO,  were  checked  by  using  the  defini- 

p 

tions  +  yc 

«  =  W/frw/p^  ,  (42) 

Y*       =         /H7W/pw^>      ,  (43) 

and    fAt)+=        1"^  ,„  W 

with  the  definitions  of  the  various  dimensionless  quantities 
to  derive  the  comparison  equations 

U+   =  -|C   (U)  -  (45) 

y+    r  ±^(/--n)   ,  (46) 

and    M*Jz   J  -£ti£  (C^)   .  (47) 

J.  /    4*  \ 

Figure  2  shows  the  comparison  of  u'  and  ~  (&t  )  for  prob- 
lem 429B  With  the  universal  velocity  and  temperature  profiles. 
The  values  of  velocity  coincide  almost  exactly  up  to  y^  =100, 
then  they  start  to  diverge  slowly  below  the  universal  profile. 
The  values  of  temperature  fall  slightly  above  the  universal 
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temperature  profile.  The  value:  -   pro;  1<  :n  397  concide  exactly 
with  those  of  problem  429B  excepi  oho  temperature  starts  to 
diverge  above  the  values  of  problem  4293  at  y"*=300.   Volume 
heat  sources  had  negligible  effect  on  the  velocity  and  temper- 
ature profiles . 

The  differences  in  the  values  of  these  two  problems  and 
problem  3  of  Dr.  Ojalvo  are  attributed  to  the  improved  expre.sr 
sion  for  the  relative  viscosity,   A  decrease  in  F  from  1000  to 
.5,  that  is  an  increase  in  the  volume  heat  sources^  had  very 
little  effect  on  the  profiles, 

All  of  the  problems  in  part  C  with  a  Prandtl  number  of  1 
or  10  and  with  a  volume  heat  source  parameter  F~10  are  plotted 
in  Figure  3.?  which  shows  the  variation  of  Nusselt  number  with 
Reynolds  number.   The  Dittus-Boelter  equation'  for  pure  forced- 
convection  turbulent  heat  transfer  is  also  plotted*  In  both 
Prandtl  number  comparison.:,  the  results  of  the  present  analysis 
compare  almost  exactly  with  the  curves  of  the  Dittus-Boelter 
equation.   Volume  heat  sources  had  negligible  effect  on  the 
variation  of  Nusselt  number  with  Reynolds  number. 

Figures  4  through  10  show  results  for  the  problems  in  tur- 
bulent flew,  combined  forced  and  free  convection.   The  Nusselt 
number  is  plotted  against  ftaylej  '  numb  :r  in  Figure  4.   The 


Nusselt  number  is  increased  less  than  an  order  of  magnitude  as 
Re'-  Is  increased  an  order  of  magnitude.   Increasing  the  Prandtl 
number  has  less  effect  on  increasing  the  Nussult  number  than 
increasing  Re*.   Increasing  the  Rayleigh  number  had  no  effect 
on  Nusselt  number.  Decreasing  P  had  a  very  little  increasing 
effect  on  Nu. 

Several  test  problems  were  run  to  determine  the  limitations 
of  the  computer  for  Re*=10^.   No  solutions  were  obtained  on  the 
IBM  1620  computer  for  Re*=lCW,  F-.l  or  less,  and  Ra=625  or  more. 
Only  by  decreasing  Ra  or  increasing  either  Re*  or  P  would  the 
computer  give  a  solution.   A  large  Ra  in  combination  with  a 
small  F  (a  large  volume  heat  source)  would  physically  appear 
to  be  opposing  each  other.   Convective  forces  predominate  in 
a  large  Ra,  whereas  with  a  small  F,  considerably  more  heat  is 
being  generated  than  is  being  convected  downstream. 

Figure  5  shows  the  pressure-drop  parameter  C  plotted 
against  Rayleigh  number.   C  increased  approximately  an  order  of 
magnitude  as  Re*  increased  an  order  of  magnitude.   For  Rayleigh 
number  less  than  625  j  Ra,  P,  and  Pr  had  no  effect  on  C.   For 
values  of  Ra  over  625  and  F}1,  increasing  Pr  decreases  C  and 
increasing  Ra  increases  C.   If  F<1,  increasing  Pr  increases  C 
and  increasing  Ra  decreases  C.   These  curves  show  that  C  changes 
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slope  for  a  value  of  P  approximately  equal  to  1.   For  a  given 
Re*  and  Pr  the  point  at  which  C  changes  slope  is  a  function  of 
Ra.  For  small  values  of  F  and  Re*,  and  large  values  of  Ra,  C 
would  appear  to  go  negative  verifying  Hallmarks  results  for 
laminar  flovH.  As  Re*  is  increased  the  pumping  action  of  large 
heat  sources  is  reduced  considerably. 

Figures  6,  Js    and  8  shew  the  variation  of  mixed- mean- to- 
wall  temperature  difference  with  Rayleigh  number.   Increasing 
Re*  an  order  of  magnitude  approximately  decreases  ty      an  order 
of  magnitude.   Increasing  Ra  has  no  effect  on  (j)  .   Increasing 
Pr  has  less  effect  on  decreasing  (j)  m  than  increasing  Re*. 
Finally  Increasing  F  has  the  most  effect  on  decreasing 


m 

As  shown  in  Figure  8,  (j)   passes  through  zero  for  F  approxi- 
mately equal  to  1.   If  (j)  m  is  plotted  against  l/F,  as  is  shown 
in  Figure  93    it  is  seen  that  a  variation  In  Ra,  Pr,  or  Re*  has 
little  effect  on  this  crossover  point.   Hallman's  curve  for 
laminar  flow,  Ra-0  Is  plotted  for  comparison. 

When  rt)  m=0,  Nu=oo,  according  to  equation  D-3.   Figure  10 
Shows  Nusselt  number  plotted  against  l/F  for  Ra=0,  Re*=lo3,  and 
Pr  of  1  and  10.   An  infinite  discontinuity  occurs  whenever 
(bm-0.   From  equation  D-3  it  is  seen  that  Nu=0  for  F=l.   Hall- 
man's  plot  for  laminar  flow,  Ra=0  is  shown  for  comparison. 
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Representative  velocity  and  temperature  profiles  are  shown 

■'  v 

in  Figures  11  through  21. 

The  effect  of  varying  Re*  is  shown  in  Figures  11  through 

16  for  Pr--10,  Ra=625,  and  for  various  values  of  F.   As  Re* 

o  b 

increases  from  10J  to  10  '",  the  velocity  in  the  center  of  the 

tube  becomes  less  positive  and  the  velocity  near  the  wall  be- 
comes more  positive.  As  Re-  is  increased,  the  temperature 
profiles  flatten  out  and  the  slope  becomes  higher  at  the  wall. 
Decreasing  F,  as  shown  in  Figure  21,  flattens  the  velocity  pro- 
file, whereas  the  temperature  profile  decreased  considerably 
and  passes  through  zero  at  F=l.   The  effect  of  increasing  F  on 
the  velocity  and  temperature  profiles  was  not  as  is  predicted 
in  appendix  F.   The  action  is  due  to  the  increased  values  of 
Re*  investigated  in  the  turbulent  flow  problems. 

Figures  17  and  18  show  the  effect  of  various  Prandtl  num- 

o        h 

bers  for  Re^'rlCK,  RarlO  ,  and  F=10.   As  Prandtl  number  increases, 

the  velocity  and  temperature  profiles  decrease  more  in  the  cen- 
ter of  the  tube  than  near  the  wall,  thus  becoming  more  flat. 
The  effect  of  a  variation  of  Ra  for  Re*=lCW,  Pr-10,  and 
Fr.5  is  shown  in  Figures  19  and  20.   As  Ra  is  increased  from  0 
to  10  there  is  very  little  decrease  in  the  velocity  and  tem- 
perature profiles  at  the  center  of  the  tube. 
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CHAPTER  V 
SUMMARY  AND  CONCLUSIONS 
The  effect  of  the  parameters  Re*,  Pr,  Ra,  and  F  on  Nu,  C, 
$m>  ^*   anc^  $  Profiles  is  summarized  in  Table  2. 

TABLE  2 

SUMMARY  OP  EPPECTS  IN  COMBINED 
FORCED  AND  FREE  CONVECTION 


Change  on  Nu 
for 


F>1 

increasing 


Ra<625      rq>625 


4>m    u 

Velocity  Temperature 

profile   difference 
profile 


Re*  increase  increase  increase  decrease  flattens  flattens 

curve     curve  and 
makes  <p 
approach 
zero 

Pr   increase  no  effect  decrease  decrease  lowers    same  as 

velocity  Re* 
in  cen- 
ter of 
tube, mak- 
ing it 
less  pos- 
itive 

Ra   no  effect  no  effect  increase  slight  increase  decreases 

decrease  in  center  in  center 

at  high  of  tube   of  tube 
Ra 
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Table  2  (Continued) 

Change  on  Nu  C  ft)rn  U  <f> 

for  Ra<625    Ra>625  Profile      Profile 

Pal 

increasing 

Re*     0      increase  increase  decrease  flattens  negligible 

curve     effect 

Pr      0     no  effect  decrease  decrease  no  effect  negligible 

effect 

Ra      0     no  effect  no  effect  no  effect  no  effect  no  effect 

F<  1 

increasing 

Re*  increase  increase  increase  decrease  flattens  same  as 

curve     F  >  1 

Pr   increase  no  effect  increase  decrease  flattens  same  as 

curve     F  >  1 

Ra  no  effect  no  effect  decrease  no  effect  flattens  negligible 

curve     effect 
slightly 

decreasing 

F    slight   no  effect  decrease  increase  flattens  increase, 
increase  curve     passing 

through 
zeroy  P=l 
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Conclu signs 

1.  All  test  problems  checked  sufficiently  close  to  previous 
known  results  to  insure  the  validity  of  the  program. 

2.  The  new  expression  used  for  the  relative  viscosity  enabled 
the  universal  velocity  and  temperature  profiles  to  be  approached 
almost  exactly  in  pure  forced-convection  turbulent  heat  transfer 

y   problems  with  negligible  volume  heat  sources. 

3.  Volume  heat  sources  had  negligible  effect  on  the  velocity 
profile,  but  had  considerable  effect  on  the  temperature  profile 
as  would  be  expected  from  the  development  of  the  equations. 

4.  The  program  was  unable  to  solve  problems  involving  high 
Rayleigh  numbers  and  volume  heat  sources  for  low  values  of  Re*. 

5.  Increasing  the  various  input  parameters  tend  to  flatten 
the  velocity  and  temperature  profiles  as  would  be  expected. 

6.  The  effect  of  volume  heat  sources  was  reduced  considerably 
by  the  influence  of  turbulent,  flow. 
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Figure  1.   Simplified  Flow  Diagram  of  the  Calculations 
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Figure  2.   Comparsion  of  Problem  429B  and  397  Results 
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Figure    3.      Comparsion   of  Turbulent  Flow,    Pure 
Forced   Convection    Results   with 
Dittus-Boelter  Equation 
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w  based  on  angular  position 

SUPERSCRII 


? 


!  t 


differentiation  once  with  respect  to  the  independent 
variable,  "ft 


per  unit;  .  .  • 
1 ' f   per  unit  volu 

Qr       Grashof  rentier ,        ^m^3MJ>L^ 

16; 

Nu  Husselt  number,   hl>/k 

Pr  Erandtl  number,      Z?/k/*t{     -    ~*/^„ 

Ra  Raylej    h  ;       '  :  <r,  -&*£?  ^ »  -   &r  Pr 

Re  Reynolds  number,   ^M^D^U   z  wm  D^ 

Re-  Friction  P.  ._.■  Dlds  number,  Du*/s  ' 
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APPENDIX  B 
THE  EDDY  DIFFUSXVITY  OP  MOMENT! 

An  analytical  study  was  made  by  Jackson^"  of  the  edrly 
diffusivlty  of  momentum  in  a  vertical  circular  tube.  The 
ratio  of  eddy  diffusivlty  to  kinematic  viscosity,  -~j~~     > 
called  relative  viscosity  appeared  In  a  reduced  form  of  the 
Navier-Stokes  equation.   An  empirical  fit  of  experimental  data 
was  employed  to  determine  the  dimensionless  velocity  gradient 
appearing  in  the  expression  for  -™  .1^11 

The  result  of  his  study  used  in  the  present  analysis  is 
the  following  set  of  expressions  for  variation  of  relative 
viscosity  with  dimensionless  radius 

£§  z  . I  Re*-*  ( 1  -y\)  -  I  for  A.  <  7)< /  ~~  £0 

e„  -  Q  R*Ssoo  ....  1  Ar  0  ^  ^}  <  -L 

V  '  (      to 

The  re  lults  compare  favorably  with  the  experimental  and 
analytical  results  of  others  and  appears  to  be  an  improvement 
over  the  ~~r     expression  used  in  recent  studies.   These  re- 
sults appear  as  equations  (29) j    (30),  and  (31)  in  the  present 

analysis. 
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APPENDIX  0 


TfiE  LONGITUDINAL  TEMPERATURE  GRADIENT 


The  heat  balance  for  the  steady-state  condition  in  oar 
problem  states  that  the  rate  at  which  heat  is  added  to  the 
fluid  is  equal  to  the  rate  at  which  heat  is  eonveeted  down- 
stream in  the  fluid.   In  symbolic  form  (see  Figure  C-l) , 


u 


/»/■ 


where  qw"  is  taken  positive  if  heat  is  being  removed.   Solving 
equation  (C~l)  for  dt/$z    gives 


The  mean  velocity  is  defined  as  follows: 


UlTi      E 


uzfrrdr 


8  I     ur d 


rr  o  */\ 


D 


from  v;hich 


/. 


m  rdr 


U&  D' 
8" 


Equation  (C-2)  can  therefore  be  written  as 


hi 


.in 


ro-  ft; 


(C-2) 


(C-3) 


(C-4) 


(C-5) 


Since  all  quantities  on  the  right-hand  side  of  equation  (C-5) 
are  constant  in  our  problem, 
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FIGURE  C-l 

Heat  Balance  for-  an  Element  of  Fluid 

Therefore,  the  temperature  gradient  in  the  direction  of 
flow  is  a  constant  independent  of  radial  location. 
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APPENDIX  D 
NUSSELT  NUMBER 

Nusselt  number  Nu  is  a  useful  expression  in  comparing 
our  problem  results.  It  is  function  of  the  temperature  grad- 
ient at  the  wall  to  a  constant  average  temperature  gradient-3. 

Considering  a  heat  balance  on  a  unit  length  of  tube: 

<?'"  rrb     a"  rru  =  o   umcrFI  rr_bz  0>-i) 

jr       if.  Lw  \m  (  a 

where  qw"  is  taken  positive  if  heat  is  being  removed-^.   Sub- 
stituting in  equation  (D-l)  the  equations 

%i"   =  h  0m 
and      (hm  i  1£JlBjxl__ 

gives     %_. _.  hD    <b»n    %'"   __   _  i         (D-2) 

Substituting  the  definition  of  P  In  equation  (D~2)  and  rear- 
ranging terms  gives 

Nu  =  hD/k  =   -f-  (l-F)  (D-3) 
From  this  equation,  as  (j)rn  approaches  zero,  Nu  will  ap- 
proach infinity;  and  when  P  =  1,  Nusselt  number  will  equal  -zero 
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APPENDIX  E 
REYNOLDS  NUMBER 

Reynold  number  Re  is  the  ratio  of  the  inertia  forces  to 
frictional  forces  and  is  a  measure  of  turbulence. 

Considering  a  dynamic  force  balance  on  an  element  of 
fluid,  the  following  expression  can  be  developed: 

&    -    &$.      ~      T/iik    n   A  (E-l) 


where  f  is  for  upward  flow  and  -  is  for  downward  flow  . 
Prom  equation  (13)  the  pressure-drop  parameter  C  is 

c  =     -^cC/f  +  fwl)  <E-2> 

3  £  yA  u  m 

and  substituting  equation  (E-l)  in  equation  (E-2)  gives 

r  _     -&%[**£' +(U- M  il  (b-3) 

w  —  a  9  ^j  1 1 ... 


3  2  M   u  ^ 
Using  equations 

rw  =  ( i/)2  (^/5c  ; 

and  the  definitions  of  Re5   Re*,   Ra,   and  P,   we  obtain 

C  -       *  (fa  )    -     _        /ft*   4m         ,  fE-4) 

Solving  for  Re  gives 

He- *    (l?*^ (E-5) 

*        8C  +    R«  <K/^F 

which  is  equation  (36)  in  the  Analysis, 
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APPENDIX  P 
VOLUME  HEAT  SOURCES 

The  term  oxu  accounts  for  the  volume  heat  sources  in 
our  problems.  Since  only  uniform  volume  heat  sources  are 
treated.,  qMI  is  independent  of  direction. 

The  nondimensional  volume  heat  source  parameter  is 

F  =    P™  Cp  Urn   R  (P-l) 

as  defined  by  Hallmark .   The  physical  significance  of  P  can 
be  described  as  the  ratio  of  the  thermal  energy  convected 
downstream,  per  unit  volume  to  the  heat  generated  in  the  fluid 
per  unit  volume.   Considering  a  heat  balance  on  a  unit  length 
of  pipe  with  constant  v/all  heat  addition 

and  introducing  equation  (F~l)  gives 


(P-3) 


P  r  1 

The  meaning  of  the  various  values  of  F  can  be  seen  from 


D  -  "' 


equation  (P-3) .   When  qw  =  0,  which  corresponds  to  an  insulated 
tube  wall j  then  P  =  1.   The  walls  should  be  hotter  than  the 
bulk  of  the  fluid  because  heat  is  convected  away  from  this 

6-'4 


region  less  rapidly  than  in  the  center 3   and  because  of  uniform 
volume  heat  source  generation.  Fluid  velocities  should  there- 
fore tend  to  be  higher  in  the  vicinity  of  the  wall  than  near 
the  center  line. 

When  F  < lj    it  means  that  heat  is  being  removed  at  the 
walls  because  more  heat  is  being  generated  internally  than  is 
convected  downstream.   The  velocity  at  the  wall  should  be  slow- 
er than  the  velocity  at  the  center  line  in  this  case  since  heat 
Is  being  removed  and  the  density  of  the  fluid  is  higher  at  the 
wall. 

When  F>1,  it  means  that  heat  is  being  added  at  the  walls 
because  heat  is  being  convected  downstream  faster  than  it  is 
being  generated.   The  velocity  in  this  case  should  rise  faster 
near  the  walls  than  near  the  center  line  since  it  is  hotter  at 
the  walls. 

F  =  0  will  occur  whenever  there  is  no  net  through-flow 
(um  =  0) 3    no  flow  (u  ~   0) j    or  constant  wall  temperature  (AzO) . 
An  investigation  is  being  conducted  on  the  general  problem 
with  a  constant  wall  temperature  by  Page1?. 

F  z  oo  would  correspond  to  no  volume  heat  sources.  An 
investigation  is  being  conducted  on  this  general  problem  by 
Griffin16. 
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